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ABSTRACT 

We present Infrared Space Observatory (ISO) observations of ^^OH and ^^OH 
toward Sgr B2 with a spectral resolution of ~35 km s~^. The OH J=5/2^3/2 
and J=3/2— i>l/2 rotational lines of the ladder are seen in emission while 

the cross-ladder transitions (from the ^113/2 J=3/2 to the J=l/2, 3/2 and 5/2 
levels of the ^IIi/s ladder), and the ^Ila/s J=5/2^3/2 and J=7/2^5/2 lines are 
detected in absorption. The ^^OH ^Hg/s J=5/2^3/2 A-doublet at ~120 /im is 
also observed in absorption. All OH A-doublets are resolved (except the ~98 
/im) and show, in addition to the strong absorption at the velocity of Sgr B2, 
several velocity components associated to the gas surrounding Sgr B2 and to the 
foreground clouds along the line of sight. 

No asymmetries in the line intensities of each doublet have been observed. 
We have modeled the observations using a non-local radiative transfer code and 
found that the OH absorption/emission must arise in a shell around Sgr B2 not 
resolved by the ISO/LWS beam. The gas density is moderate, with upper limits 
of 10^ cm~^ and ~300 K in temperature. The OH abundance is high, (2-5) x 10^^. 
We argue that a widespread photon dominated region explains the enhancement 
of OH abundance. 

Subject headings: infrared: ISM: lines and bands — ISM: individual (Sgr B2) 
— ISM: molecules — line: identification — molecular data — radiative transfer 
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cially the PI countries: France, Germany, the Netherlands and the United Kingdom) and with participation 
of ISAS and NASA. 
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1. Introduction 

The Sagittarius B2 complex represents an interesting burst of massive star formation in 
the Galactic Center (GC) and may be representative of other active galactic nuclei. Large 
scale continuum maps from radio to far-IR wavelengths show that Sgr B2 is the brightest 
and among the most massive clouds of the GC (Cox & Laureijs 1989, Pierce-Price et al. 
2000; Scoville, Solomon & Penzias 1975). The Sgr B2 M source is the brightest far-IR 
condensation of the complex with a diameter of ~40" at 100 //m (Goldsmith et al. 1992) 
and has also the largest gas column density (Lis & Goldsmith 1989). These studies have 
shown that the core is embedded in an extended and clumpy cloud of warm gas. However, the 
observed rich chemistry in the gas in front of Sgr B2 and its possible heating mechanisms are 
far from settled. Low velocity shocks have been invoked to explain the enhanced abundances 
of some species such as NH2 or NH3 (Flower, Pineau des Forets & Walmsley 1995) which 
are not observed in more quiescent regions. On the other hand, UV radiation can have an 
important effect on the gas in the outer layers. This radiation may be provided by evolved 
stars and by young massive stars in the envelope of Sgr B2 itself (Martm-Pintado et al. 
1999). In addition, the emission of several ions with high excitation potential (such as [OIII] 
and [NIII]) is extended in the Sgr B2 region (Goicoechea et al. 2003, in preparation). Thus, a 
widespread ionized component producing photodissociation regions (PDRs) in the envelope 
is also possible. In any of these scenarios or in a combination of both, 0-bearing species 
such as H2O, OH, H3O+ and atomic oxygen are decisive for the thermal balance of the 
dense molecular gas (Neufeld et al. 1995). In particular, the hydroxyl radical, OH, has been 
predicted to be abundant in both scenarios. 

The structure of the OH levels is depicted in Figure 1. The fundamental rotational 
A-doublet (~119 /xm) was first detected in the direction of Sgr B2 in absorption against the 
thermal dust emission by Storey, Watson & Townes (1981). Each line of the doublet is so 
optically thick that absorbs completely the continuum radiation avoiding any check with the 
models. These lines are detected over a large area (9'x27') around Sgr B2 (Cernicharo et al. 
1999). A detailed study of the far-IR spectrum of OH has been carried out mainly in the 
Orion shocked region (Watson et al. 1985; Viscuso et al. 1985; Melnick et al. 1987, 1990; 
Betz & Boreiko 1989) but little is known toward GC clouds such as the gas in the line of 
sight toward Sgr B2. 

Prom the theoretical point of view. Offer & van Dishoeck (1992; hereafter OfD) and 
Offer, van Hemert & van Dishoeck (1994; hereafter OfHD) have made detailed calculations 
of the collisional cross sections between OH and H2 and of the expected intensity of the OH 
far-IR lines. They predicted strong asymmetries in the intensity of the OH A-doublets due 
to important asymmetries in the collisional rates between OH and para-H2. 
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In this letter, we present high-resolution observations of several OH A-doublets toward 
Sgr B2 involving levels up to ~420 K. All lines inside the ladder are seen in emission 

while all the remaining lines are detected in absorption. The data have been modeled using 
a non-local radiative transfer code in order to constrain some important physical parameters 
of the warm outer layers of Sgr B2 and of the foreground gas. 

2. Observations, Data Reduction and Results 

Most of the pure rotational lines of OH are in the far-IR wavelength coverage of the 
Long- Wavelength Spectrometer (LWS; Clegg et al. 1996) on board ISO (Kessler et al. 
1996). The ISO data base observations used in this Letter are 32201429, 32600502, 46201118 
46201123, 46701803, 47600809, 47608099, 49401705, 50400823, 50601013, 50600603, 50700610, 
and 84500102. The LWS/Fabry-Perot (FP) spectral resolution is ~0.015 /xm (~35 km s"^) 
and allows resolving the A-doubling of OH and its isotopes. Several OH rotational hnes 
(intra and cross-ladder) were searched toward Sgr B2 M. Data are shown in Figure 1. The 
OH J= 5/2<— 3/2 cross-ladder doublet at ~34 //m was observed with the Short- Wavelength 
Spectrometer (SWS; de Graauw et al. 1996) with a spectral resolution of 1500 (ISO observa- 
tion 28702002). One component was observed with the SWS/FP with a velocity resolution 
of ~10 km s~^ (ISO observation 46001217). When compared to the grating observation, the 
mean LWS/FP continuum flux of each line deviates by at least 20%, similar to the grat- 
ing flux calibration errors (Swinyard et al. 1998). The systematic uncertainties are always 
larger than the statistical errors associated with the line signal-to-noise ratios (S/N). The 
best calibration is obtained by correcting the continuum level of the FP spectrum in order to 
coincide with the grating level. A polynomial baseline was fitted to each spectra and adopted 
as the continuum level. The effective SWS/FP aperture in the ~34 //m range is 17"x40", 
while the LWS aperture is circular and ~80" in diameter. All the SWS and LWS products 
were processed with version 10.1 of the pipeline and analyzed using the ISO spectrometers 
data reduction package ISAP.^ 

Sgr B2 M has been the target of several spectroscopic observations in the far-IR. Molec- 
ular features have been always found in absorption. However. OH shows a different behavior. 
All rotational doublets connecting the ^^OH ^Il3/2 ground-state with the rotational levels of 
the ladder (~34, 53, 79 and 119 fim) arc observed in absorption while the ^®0H ^Hi/2 

intra-ladder rotational transitions at ~98 and ~163 /xm are observed in emission. 



■^The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instruments 
Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC,MPE, RAL, and SRON. 
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The ~79 and ~53 /xm A-doublets have maximum absorption depths at the velocity 
of Sgr B2 M (VL5ij~+50 km s""^). Nevertheless, there is still a contribution from more 
negative velocity components associated with the "spiral arms clouds" in the line of sight 
and with the gas surrounding the GC, which we divide in two components {c^ —85 and 
km s~^). The n(H2) range found in these clouds varies from 10^ to 10^ cm~^ (Greaves 1995). 
The absorption depths at Vl5r— -85 km s^^ are ~30% and ~15% for the ~79 and ~53 /xm 
doublets respectively. These components are thought to arise in the molecular gas within 1 
kpc of the GC (Scoville 1972). Assuming that the absorption at these velocities is optically 
thin and considering only radiative excitation (i.e., low Te^;), we derive x{OH) ~ 10~^ for 
A^(//2)^(5-10)xl02i cm-2 (see Neufeld et al. 2000 for a distribution of N{H2) in the line 
of sight). The opacity of the OH lines at Vlsr—^ km s~^ allows also to derive x{OH) ~ 
5 X 10^^ for A^(iJ2)— 2.5x10^^ cm~^. This component is believed to arise in the diffuse gas 
at galactocentric radii 3-8 kpc and in the Solar vicinity (Greaves & Williams 1994). These 
values are higher than the typical OH abundances of 10~^ found in diffuse and translucent 
clouds (Van Dishoeck & Black 1986). 

The ~84 fim doublet and the 34.603 /xm line component (in absorption) and the ^Hi/2 
doublets at ~98 and ~163 fim (in emission) are strictly centered at the velocity of Sgr B2 
M and are decisive to understand the excitation scheme of the molecule. The ISO data 
also provide upper limits to the emission/absorption of the J=3/2— >^5/2 (~96 fim) and 
J=5/2— >5/2 (~48 /xm) cross-ladder transitions and to the ^n3/2 J=9/2^7/2 (~65 /xm) and 
^ni/2 J=7/2— >5/2 (~71 iim) A-doublets of 1.5% of the continuum emission. Finally, we 
present the detection of the fundamental ^^OH ^n3/2 J= 5/2-^3/2 hnes at ~120 /xm. Apart 
from the feature associated with Sgr B2, there is another absorption at Vlsr —0 km s~^. The 
derived ^*^0H column density is (2.5±0.5) x 10^^ cm~^, implying an ^^OH/^®OH abundance 
ratio in this component of ~500, similar to the terrestrial value. 



3. Modeling Sgr B2 and Discussion 

The typical OH abundance in dense quiescent clouds is (0.1-1) x 10~^ with OH/H2O 
ratios in the range 1 to 10~^ (Bergin et al. 1995). According to the chemical models, the 
major contribution to the enhanced OH abundance comes from regions where water vapor is 
being rapidly photodissociated (Sternberg & Dalgarno 1995) or regions in which shock waves 
play a role (Draine et al. 1983). In the case of Sgr B2, the innermost regions of the cloud 
are completely hidden in the mid and far-infrared due to the large dust opacity. Almost all 
the observed far-IR OH come from the external layers of the Sgr B2 cloud, i.e., from the 
surrounding gas at high Tk and from the cold dark clouds and diffuse gas in the hne of sight. 
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In order to estimate the OH abundance and the physical conditions leading to the 
observed OH line emission in the Vlsb. — +50 km component, we have modeled the first 
20 OH rotational levels using the code developed by Gonzalez- Alfonso & Cernicharo (1993). 
The hyperfinc structure of OH is not included in the model. The collisional cross sections 
of OfHD have been used (kindly provided in electronic form by E.F van Dishoeck). We 
have adopted a spherical geometry for a cloud consisting of two components: a uniform core 
with a diameter of 25" (3.2x10^^ cm for an assumed distance to Sgr B2 of 8.5 kpc), with a 
dust temperature of 30 K (see Goicoechea & Gernicharo, 2001) and opacity at 80 //m of 5 
(ta = Tso * (80/A(/xm)), and a shell of variable thickness and distance to the core. The shell 
was divided into 14 layers. All molecular transitions in the core are thermalized to the dust 
temperature due to the large dust opacity in the far-IR. In order to check the sensitivity of 
the results on the physical parameters, n{H2), and N{OH) were varied from 10^ to 10^ 
cm~^, 40 to 600 K, and 1x10^^ up to 1x10^^ cm~^ respectively. 

The observations put some constraints on the size and density of the shell. Excitation 
temperatures for the cross-ladder and ^11^/2 transitions have to be below the dust temperature 
in order to reproduce the observed line absorptions. In addition, the shell thickness and 
its distance to the core cannot be large compared to the core size, otherwise limb effects 
introduce important re-emission that could cancel the absorption produced by the gas in 
front of the core. The results for some models are reproduced in Figure 2. We have found 
that if the shell is placed at some distance from the core, then it is difficult to reproduce the 
observations as almost all lines appear in emission within the LWS beam (80") . Even if the 
shell is contiguous to the core but large in thickness, lines also appear in emission. Models 
in Fig. 2 correspond to a total size of 42" (5.3x10^^ cm) for the core-|-shell cloud. 

The model labelled Mi corresponds to Tx=40 K, n(//2)=10^ cm-^ and x{OH) = 2 x 
10""^. The cross-ladder and ^11^/2 transitions are in absorption while the '^Il.\/2 intra-ladder 
transitions are predicted in emission. However, although the intensity of the absorbing lines 
is found to fit reasonably well the observations, emission lines are underestimated. If the 
OH abundance is increased in order to reproduce the intensity of the emitting lines then, 
absorbing lines are poorly fitted as limb effects start to dominate. This effect is already seen 
in the J=l/2<— 3/2 cross-ladder transition at ~79 /xm. In particular, the ^Hi/2 J=5/2^3/2 
transition at ~98 /xm is strongly underestimated by model Mi. If the gas temperature 
increases, collisions start to pump the OH levels and the intensity of the ~98 //m line is 
better fitted. However, a new problem arises as strong asymmetries in the line intensities 
of each A-doubling line do appear at high temperature (no asymmetries have been observed 
within the S/N ratio of the data ). These asymmetries were already predicted by OfD and 
are due to a strong parity change propensity rule introduced by collisions with para-H2. 
They can only be suppressed if the H2 ortho-to-para (OTP) ratio is higher or if radiative 
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pumping dominates. We have adopted a OTP ratio of 3. Model M2 in Fig. 2 corresponds 
to Tj^=150 K with otherwise the same parameters as Mi. The intensity of the ~98 /xm 
doublet is considerably enhanced with respect to the low temperature model. The intensity 
asymmetry in this doublet is not visible in the data as the doublet is unresolved with the 
LWS/FP spectrometer. The effect of the IR pumping strongly reduces the asymmetries that 
could be present without the optically thick core. Finally, model M3 (thick lines in Figure 
2) represents Tk=SOO K, n{H2)=5xlO^ cm~^ and x{OH) = 3 x 10~^. Again, the strong 
asymmetry of the ~98 jjim doublet is not observable due to the lack of spectral resolution. 
However, the observed total intensity for the doublet agrees well with that predicted from 
the model. For the other lines, M3 reproduces qualitatively well the observed OH pattern, 
except in the line wings in the 119 and 79 jim. doublets. Nevertheless, cold foreground gas 
(not included in the model) absorbs (like in Mi) at these wavelengths, reducing the emission 
wing effects in these lines. 

From our models it is clear that the gas surrounding Sgr B2 has to be warm (see 
Ceccarelli et al. 2002) and must have a moderate H2 density. The OH abundance is rather 
high, (2-5) X 10^^. Assuming that M3 represents a reasonable fit to the OH data, we derive 
N{}^OH) = {<6±2)xW^ cm-2 for Sgr B2. Hence, the ^^O/^^O isotopic ratio is 240-280, in 
excellent agreement with the value obtained by Bujarrabal, Cernicharo & Guelin (1983) from 
the OH 18 cm lines and a factor ~2 lower than in the V lsr —0 km s~^ component. Higher 
Tk could even be possible but the actual coUisional rates introduce very strong intensity 
asymmetries in the cross-ladder transitions which are not observed. If the gas temperature is 
higher, then its density has to be lower than the value used in M3 (5x10^ cm~^). The reason 
is that higher temperatures produce strong emission in the OH high excitation transitions 
(~71 and ~65 jim. for example, see Figure 2) which is not observed. For Ti^=600 K, n{H2) 
has to be decreased to (1-2) x 10^ cm^^ and x{OH) increased to (0.5-1) x 10^^ to obtain a 
crude fit the ISO data. A temperature gradient across the shell, going from T/^=40 K in the 
innermost regions to T;^=600 K in the external layers is probably a better representation to 
the physical structure of the edge of the Sgr B2 cloud. 

The inferred ^^OH column density in the shell, (1.5-2.5) x 10^^ cm'^, and that of H2^0 
determined also from LWS/FP data (Cernicharo et al. 1997; 2002, in preparation) leads 
to an OH/H2O=0.1-l abundance ratio. For comparison, shock models of Flower et al. 
predict ~10~^, while models of dense dark cores predict ~10~^. On the other hand, the 
maximun OH abundance in PDR models is ~10~^ with OH/H2O ratios close to 10. The 
large abundance found in Sgr B2 suggests that its external shells are illuminated by a strong 
UV field producing a PDR, although low velocity shocks (v^ ~30 km s~^), unresolved by 
the ISO data, may be also present. Many species formed during the evolution of the cold 
gas in Sgr B2 are being now reprocessed in these regions. This warm gas is poorly traced 
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by millimeter and submillimeter observations, but it represents the strongest contribution 
to the absorption/emission features in the far-IR spectrum of Sgr B2. In addition to ^^OH 
and ^^OH, we have searched for several related species such as 0H+ or H2O"'". No lines from 
these species have been found. 

The future far-IR instruments on board next generation telescopes such as the Herschel 
Space Observatory will allow a fast mapping of many OH rotational lines in galactic and 
extragalactic sources making OH a very useful tool in deriving the physical properties of GC 
giant molecular clouds. This could help in the analysis of starburst galaxies like Arp 220. 

We thank Spanish DGES and PNIE for funding support under grant ESP98-1351-E and 
PANAYA2000-1784 and CM. Walmsley, V. Bujarrabal and J. Fischer for their comments 
on the manuscript. We also thank our referee, E.F. van Dishoeck, for her useful comments 
and suggestions. JRG acknowledges UAM for a pre-doctoral fellowship. 
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Fig. 1. — ISO observations of ^^OH and ^^OH toward Sgr B2 M and rotational energy 
diagram of OH. The two ladders, ^113/2 and produced by the spin orbit interaction are 

shown. The A-doubhng splitting of each rotational level has been enhanced for clarity. The 
hyperfine structure of the rotational levels is not shown. The observed rotational transitions 
are indicated by arrows and the observed ISO line profiles are shown in boxes. The intensity 
scale corresponds to the continuum normalized flux and the abscissa to the wavelength in 
microns. When lines are centered only at Sgr B2 velocities (~+50 km s~^), the LSR velocity 
scale is not shown. All A-doublets are LWS/FP data excepting the 34.6 /xm doublet for 
which SWS grating data (A/AA ~1500) are shown. One component has been observed with 
the SWS/FP (A/AA ~30000) and data are shown in the inset of the corresponding box. The 
^^OH ^113/2 J=5/2<— 3/2 doublet is shown in the upper inset. 
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Fig. 2. — Results from selected models of Sgr B2 component (Vlsr — +50 km s^^) dis- 
cussed in the text. The central panels correspond to the cross-ladder transitions of OH. 
Line intensities have been normalized to the continuum dust emission. The wavelength 
of the transitions is indicated at the right of each panel. The model labelled Mi corre- 
sponds to Tfc=40 K, n(H2)=10^ cm~^, xoif=2xl0~^. The model labelled M2 corresponds to 
Tfc=150 K, n(H2)=10^ cm-^ xoff=2xl0-^ The model labelled M3 corresponds to Tfc=300 
K, n(H2)=5xlO^ cm"^, xoi?=3xl0~^. The velocity resolution in the modeled hne profiles 
is 1 km s~^. 
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